Comparing maps

library(Racmacs)

Different maps can be compared using the procrustesMap() function where one map is compared to another.

The example map

We’ll read in the H3N2 2004 map as an example map

map <- read.acmap(system.file("extdata/h3map2004.ace", package = "Racmacs"))
view(map)

Comparing optimization runs

We can use the procrustesMap() function to compare the top 3 optimization runs.

pc_run12 <- procrustesMap(
  map                            = map,
  comparison_map                 = map,
  optimization_number            = 1,
  comparison_optimization_number = 2
)

view(pc_run12)
pc_run13 <- procrustesMap(
  map                            = map,
  comparison_map                 = map,
  optimization_number            = 1,
  comparison_optimization_number = 3
)

view(pc_run13)